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The  paper  studied  the  performances  of  parallel  thermoelectric  generator  (TEG)  by  theoretical  analysis  and 
experimental  test.  An  analytical  model  of  parallel  TEG  was  developed  by  theoretical  analysis  and  calcu¬ 
lation,  based  on  thermodynamics  theory,  semiconductor  thermoelectric  theory  and  law  of  conservation 
of  energy.  Approximate  expressions  of  output  power  and  current  of  parallel  TEG  were  deduced  by  the 
analytical  model.  An  experimental  system  was  built  to  verify  the  model.  The  results  indicate  that  only 
when  all  of  the  thermoelectric  modules  (TE  modules)  in  the  parallel  TEG  have  the  same  inherent  param¬ 
eters  and  working  conditions,  the  parallel  properties  of  the  TEG  are  the  same  as  that  of  common  DC 
power.  The  existence  of  contact  resistance  is  just  like  the  increase  of  the  TE  module’s  internal  resistance, 
which  leads  to  the  deceases  of  output  power.  The  thermal  contact  resistance  reduces  the  output  power  by 
reducing  the  temperature  difference  between  the  two  sides  of  the  thermocouples.  The  results  derived 
from  the  model  are  basically  consistent  with  the  experimental  results,  the  model  is  suitable  for  the  per¬ 
formance  researching  and  designing  of  parallel  TEG. 
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1.  Introduction 

A  TEG  produces  a  voltage  when  there  is  a  temperature  differ¬ 
ence  between  the  hot-side  and  cold-side  of  TEG  by  thermoelectric 
effect,  which  means  that  the  thermal  energy  is  converted  directly 
into  electric  energy  by  thermoelectric  effect.  Thermoelectric  effect 
includes  Seebeck  effect,  Peltier  effect  and  Thomson  effect,  it  also 
accompanies  with  other  effects,  such  as  Joulean  effect  and  Fourier 
effect. 

Thermoelectric  generation  is  a  technology  for  directly  convert¬ 
ing  thermal  energy  into  electrical  energy,  it  has  no  moving  parts, 
is  compact,  quiet,  highly  reliable  and  environmentally  friendly. 
Because  of  these  merits,  it  is  presently  becoming  a  noticeable 
research  direction  [1]. 

In  the  aspect  of  system  analysis  and  optimization,  Rowe  and 
Gao  [2]  provided  a  practical  procedure  for  optimizing  TE  module 
geometry  guided  by  the  “economic  factor”.  Chen  et  al.  [3]  explored 
the  performance  of  the  TEG  with  multi-element  thermoelectric 
equipment  by  assuming  that  the  heat  transfer  between  the  ther¬ 
moelectric  generator  and  the  external  heat  reservoirs  obeys  the 
linear  phenomenological  law  used  in  the  irreversible  thermody¬ 
namics.  Palacios  et  al.  [4]  showed  a  methodology  to  extract 
thermoelectric  internal  parameters  from  the  information  provided 
by  performance  curves,  allowing  scientists  to  predict  the  perfor¬ 
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mance  of  the  module  at  any  working  condition.  Chen  et  al.  [5]  built 
a  model  of  a  two-stage  semiconductor  thermoelectric-generator 
with  external  heat-transfer.  Yamashita  [6,7]  built  thermal  rate 
equations  by  taking  the  linear  and  non-linear  components  in  the 
temperature  dependences  of  the  Seebeck  coefficient,  the  electrical 
resistivity  and  thermal  conductivity.  Gao  and  Rowe  [8]  fabricated  a 
novel  tube-shape  TE  module  from  four  ring-shaped  thermo  ele¬ 
ments,  and  established  an  analytical  model  of  the  tube-shape  TE 
module,  the  result  indicated  that  a  tube-shape  TE  module  could 
achieve  similar  performance  to  that  of  a  conventional  plate-like 
module,  and  had  an  advantage  in  applications  when  heat  flowed 
in  a  radial  direction. 

In  the  aspect  of  thermoelectric  technology  application  in  indus¬ 
try,  Hsu  et  al.  [9],  Masahide  et  al.  [10],  Yodovard  [11],  Yu  and  Chau 
[12],  and  Xu  et  al.  [13]  applied  the  thermoelectric  generation  tech¬ 
nology  to  the  environment-friendly  vehicle,  both  waste  heats  from 
vehicle  exhaust  and  engine  were  recovered  into  electrical  power  to 
provide  auxiliary  power  for  the  vehicle.  Chen  et  al.  [14]  presented  a 
system  efficiency  analysis  and  impact  evaluation  of  the  application 
of  thermoelectric  power  cycle  to  thermal  energy  systems,  espe¬ 
cially  Combined  Heat  and  Power  production.  In  order  to  investigate 
viability  and  further  performance  of  the  TEG  for  waste  heat  recov¬ 
ery  in  industry  area,  Gou  et  al.  [1]  and  Niu  et  al.  [15]  constructed  a 
low-temperature  waste  heat  thermoelectric  generator  setup,  the 
testing  results  and  discussion  showed  the  promising  potential  of 
using  TEG  for  low-temperature  waste  heat  recovery,  especially  in 
industrial  fields. 

In  the  aspect  of  new  thermoelectric  materials  research  and 
development,  nowadays,  a  large  number  of  works  concerning  ther- 
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moelectricity  focus  on  how  to  improve  heat-to-electricity  conver¬ 
sion  efficiency  of  thermoelectric  materials,  that  is,  how  to  improve 
the  thermoelectric  figure  of  merit  ZT  of  these  materials.  Nanotech¬ 
nology  [16,17],  novel  technology  in  solid  state  physics  and  semicon¬ 
ductor  physics  [18,19]  are  employed  to  explore  this  exciting  field. 

Despite  these  promising  results,  wide  application  of  TEG  has  yet 
been  limited  because  of  its  relatively  low  heat-to-electricity  con¬ 
version  efficiency.  However,  compared  with  conventional  methods, 
thermoelectric  generation  technology  has  great  advantage  in  some 
individual  areas.  For  example,  in  the  area  of  low-grade  waste  heat 
recovery,  as  it  cannot  be  recycled  effectively  by  conventional  meth¬ 
ods,  TEG  appears  to  have  advantages  due  to  its  entire  solid-state  en¬ 
ergy  conversion  mode.  Furthermore,  considering  waste  heats  are 
low-cost  and  even  no-cost  resources,  the  low  heat-to-electricity 
conversion  efficiency  is  no  longer  a  major  problem,  the  focus  of 
our  concern  is  how  to  get  larger  output  power  of  a  TEG. 

As  the  output  power  of  the  TEG  composed  of  only  one  TE  module 
is  very  small,  to  obtain  higher  output  power,  it  is  feasible  to  make  a 
TEG  by  multiple  TE  modules  in  series  connection  or  parallel  connec¬ 
tion  (for  short,  named  series  TEG  and  parallel  TEG  in  this  paper, 
respectively).  Zhang  et  al.  [20]  studied  the  performances  of  ser¬ 
ies-parallel  connection  for  TE  modules  by  experiment  and  numer¬ 
ical  calculation,  the  results  showed  that  when  all  the  modules  had 
the  same  parameters,  and  worked  in  the  same  condition,  the  series- 
parallel  properties  of  the  TEG  were  the  same  as  the  series-parallel 
properties  of  common  DC  power.  However,  to  a  multiple  modules 
TEG,  it  is  hard  to  guarantee  every  TE  module  works  in  the  same  con¬ 
ditions.  For  example  when  the  thermoelectric  technology  is  used  in 
the  waste  heat  recovery  of  industry  furnaces,  due  to  the  different 
surface  temperature  of  every  part  of  furnaces,  it  will  cause  the  dif¬ 
ferent  hot-side  temperature  of  every  TE  module.  Therefore,  there 
are  more  components  in  the  multiple  modules  TEG,  and  its  perfor¬ 
mance  analysis  is  more  complex  than  the  one  module  TEG.  As  pre¬ 
vious  researches  about  multiple  modules  TEG  focused  on  the 
performance  analysis  of  the  TEG  that  all  TE  modules  are  in  the  same 
working  conditions,  when  each  TE  module  of  multiple  modules  TEG 
is  in  different  working  conditions,  it  is  necessary  to  build  an  analyt¬ 
ical  model  of  series  TEG  and  an  analytical  model  of  parallel  TEG. 

The  output  performance  of  series  TEG  has  been  studied  by  us  re¬ 
cently  [21  ].  In  this  research,  the  output  performance  of  parallel  TEG 
is  studied,  an  analytical  model  of  parallel  TEG  is  proposed  based  on 
thermodynamic  theory,  semiconductor  thermoelectric  theory,  and 
law  of  conservation  of  energy.  Approximate  expressions  of  output 
power  and  current  of  parallel  TEG  are  deduced  by  the  analytical 
model.  Compared  with  the  previous  researches,  the  model  takes 
into  account  much  more  influencing  factors,  such  as  thermal  con¬ 
tact  resistance,  contact  resistance,  and  temperature  of  heat  source 
and  heat  sink  of  every  TE  module.  Just  as  the  previous  researches, 
the  model  can  be  used  to  analyze  the  output  performance  of  paral¬ 
lel  TEG  when  all  TE  modules  work  in  the  same  conditions.  But  the 
model  can  also  be  used  to  analyze  the  output  performance  of  paral¬ 
lel  TEG  when  all  TE  modules  work  in  different  conditions,  which 
cannot  be  analyzed  by  previous  researches.  The  model  can  be  used 
for  design  and  analysis  of  parallel  TEG  accurately. 

The  purpose  of  this  paper  is  to  build  up  an  analytical  model  of 
parallel  TEG  by  theoretical  analysis,  to  investigate  how  contact  ef¬ 
fect  reduces  the  output  power  of  parallel  TEG,  and  further  to  verify 
the  model  by  experimental  test. 


2.  Model  description 

2.1.  Model  structure 

The  structure  of  a  typical  one  TE  module  TEG  is  shown  schemat¬ 
ically  in  Fig.  1,  it  consists  of  a  TE  module,  a  heat  source  and  a  heat 


sink.  ATE  module  consists  of  over  100  P-N  thermocouples,  conduc¬ 
tive  tabs,  and  two  ceramic  substrates.  The  thermocouples  are  con¬ 
figured  so  that  they  are  connected  electrically  in  series,  but 
thermally  in  parallel.  Ceramic  substrates  provide  the  platform  for 
the  thermocouples  and  the  small  conductive  tabs  that  connect 
them.  Heat  source,  thermocouples,  conductive  tabs,  ceramic  sub¬ 
strates,  and  heat  sink  thus  form  a  layered  configuration. 

When  the  temperature  of  the  heat  source  and  the  heat  sink 
cause  the  temperature  difference  between  the  hot-side  and  cold- 
side  of  thermocouples,  five  thermoelectric  effects  will  occur. 
Therefore,  a  current  I  will  flow  through  an  external  load  resistance 
rL  because  of  the  temperature  difference. 

The  structure  of  a  parallel  TEG  is  shown  schematically  in  Fig.  2. 
Compared  to  the  one  TE  module  TEG,  it  consists  of  multiple  TE 
modules,  heat  sources  and  heat  sinks,  and  only  one  load  resistance. 

The  parallel  TEG  consists  of  n  TE  modules  which  are  consecu¬ 
tively  numbered  from  1  to  n,  Tih  and  Tic  are  the  heat  source  temper¬ 
ature  and  heat  sink  temperature  of  the  ith  (1  ^  i  <  n)  TE  module, 
respectively,  Tn  and  Ti2  are  the  thermocouple’s  hot-side  tempera¬ 
ture  and  cold-side  temperature  of  the  ith  TE  module,  respectively. 

2.2.  Mathematical  modeling 

The  following  assumptions  are  made  to  simplify  the  complex 
problem  of  modeling: 

(a)  All  the  TE  modules  of  the  parallel  TEG  are  the  same  model, 
which  means  that  they  have  the  same  inherent  parameters. 

(b)  The  heat  conduction  between  thermocouples  is  ignored,  as 
the  heat  conduction  from  hot-side  to  cold-side  within  the 
thermocouples  as  axial  heat  conduction  will  be  dominant. 

(c)  The  Seebeck  coefficient,  thermal  conductivity  and  resistivity 
of  semiconductor  thermocouples  are  constants. 

(d)  On  the  assumption  that  Seebeck  coefficient  remains 
unchanged,  the  influence  of  Thomson  effect  is  not  considered. 

There  are  n  TE  modules  in  the  parallel  TEG  shown  in  Fig.  2.  To  a 
TE  module,  let  m  be  the  number  of  thermocouples,  a  (VI<-1)  be  the 
Seebeck  coefficient  of  a  thermocouple,  I<  (WI<-1)  be  the  thermal 
conductance  of  all  thermocouples,  r  (Q)  be  the  internal  resistance 
of  a  TE  module,  re  (Q)  be  the  contact  resistance  of  a  TE  module 
which  including  the  contact  resistance  between  thermocouples 
and  conductive  tabs,  the  internal  resistance  of  conductive  tabs. 
As  m,  a,  I<,  r  and  re  are  all  inherent  parameters  of  TE  modules, 
according  to  the  previous  assumption,  they  have  the  same  value 
to  each  TE  module. 
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Fig.  2.  Schematic  diagram  of  a  parallel  TEG. 


The  diagram  of  thermal  conductance  and  thermal  resistance  of 
every  part  of  a  TE  module  is  shown  schematically  in  Fig.  3. 

In  Fig.  3,  I<i h  (WK-1)  and  I<ic  (WK1)  are  the  total  hot-side  and 
cold-side  thermal  conductance  of  the  ith  TE  module,  respectively. 
Ki  (WK-1)  is  the  sum  of  thermal  conductance  of  hot  side  or  cold 
side  conduct  taps,  I<2  (WK-1)  is  the  thermal  conductance  of  hot  side 
or  cold  side  ceramic  substrate.  Pi  (KW-1)  is  the  sum  of  the  thermal 
resistance  between  the  thermocouples  and  conduct  taps  on  either 
side.  R2  (KW-1)  is  the  sum  of  the  thermal  resistance  between  the 
conduct  taps  and  ceramic  substrate  on  either  side.  Rih  (KW-1)  is 
the  thermal  contact  resistance  between  the  heat  source  and  hot- 
side  ceramic  substrate,  it  reflects  the  contact  condition  between 
the  heat  source  and  the  TE  module.  Ric  (KW-1)  is  the  thermal  con¬ 
tact  resistance  between  the  heat  sink  and  cold-side  ceramic  sub¬ 
strate,  it  reflects  the  contact  condition  between  the  heat  sink  and 
the  TE  module.  Ri3  (KW-1)  is  the  thermal  resistance  of  the  heat 


Fig.  3.  Schematic  diagram  of  thermal  conductance  and  thermal  resistance 
part  of  a  TE  module. 


Kic 


Kih 


of  every 


sink,  it  is  related  to  the  cooling  method  adopted,  the  heat  sink 
structure,  application  environment,  and  so  on. 

Since  Klt  K2,  Pi  and  R2  are  the  inherent  parameters  of  TE  mod¬ 
ule,  they  have  the  same  value  to  every  TE  module.  However,  Rih 
and  Ric  are  not  the  inherent  parameters  of  TE  module,  they  reflect 
the  contact  condition  between  the  TE  module  and  the  heat  source 
or  heat  sink,  respectively,  and  it  cannot  be  guaranteed  that  all  con¬ 
tact  surfaces  have  the  same  contact  condition  for  different  TE  mod¬ 
ules.  For  example,  when  TEG  is  used  to  convert  the  surface  heat  of 
boiler  into  electric  energy,  the  boiler  surface  will  be  the  heat  source 
of  TEG,  as  the  boiler  surface  is  not  smooth,  the  contact  condition 
between  boiler  surface  and  different  TE  module  is  not  the  same. 
So  it  is  possible  that  there  are  different  Plh  and  Ric  in  different  TE 
modules.  Besides,  as  Pl3  is  related  to  the  application  environment, 
it  is  also  possible  that  there  are  different  R&  in  different  TE 
modules. 

It  can  be  seen  form  Fig.  3  that  I<ih  and  I<ic  consist  of  many  parts  of 
thermal  conductance  and  thermal  resistance,  so  I<ih  and  I<ic  can  be 
expressed  as 

ii~  —  ir~  +  it~  +  +  ^2  +  (1) 

Ajh  A 1  A  2 

11  =  11~  +  11~  +  +  i^2  +  Ri3  +  Ric-  (2) 

AjC  Al  A  2 

For  commercial  TE  module,  the  value  of  Ki,  R i  and  R2  is  gener¬ 
ally  very  small,  the  influencing  factors  of  thermal  resistance  of  hot- 
side  and  cold-side  are  mainly  I<2f  Ri3,  Rih  and  Ric. 

It  is  assumed  that  every  TE  module  in  the  parallel  TEG  works  on 
different  conditions,  such  as  different  heat  source  temperature,  and 
different  heat  sink  temperature.  Let  Qjh  (W)  present  the  heat  that 
absorbed  from  its  heat  source  by  the  ith  TE  module  per  unit  time, 
and  Qic  (W)  present  the  heat  that  released  to  its  heat  sink  per  unit 
time. 

According  to  Fourier’s  law,  Q/h  and  QjC  can  be  expressed  as 


Qih  =  Kfli(Tfl,  —  Tn), 

(3) 

Qic  =  Kk(Ta  -  T,c). 

(4) 

In  addition,  Q,h  consists  of  three  parts  of  heat,  Pletier  heat  on  the 
hot-side  of  thermocouples,  Joule  heat  generated  by  internal  resis¬ 
tance  of  the  TE  module  as  the  current  flow  through  the  TEG  circuit, 
and  conduction  heat  from  the  hot-side  to  the  cold-side  of  the  ther¬ 
mocouples  [22,23].  So  Qih  can  also  be  expressed  as 

Q,h  =  mocTn  /,  -  (r  +  re)  +  I<(Tn  -  Ta ) ,  (5) 

where  /*  is  the  current  flowing  in  the  ith  TE  module. 

Similarly,  QiC  can  also  be  expressed  as 

Qic  =  mocTal,  +  1/f  (r  +  re)  +  K(Tn  -  Ta ) .  (6) 

Let  the  current  flowing  through  the  load  be  J,  according  to 
Kirchhoffs  current  law,  in  a  parallel  circuit  the  total  current  is 
equal  to  the  sum  of  the  currents  in  each  individual  branch,  the 
current  /  can  be  expressed  as 

/  =  /i+/2  +  . ../i  +  . ..  +  /„  =  £/,,  (7) 

i=  1 

Let  the  load  resistance  of  the  parallel  TEG  be  rL  (Q),  the  output 
power  of  the  parallel  TEG  be  P  (W),  and  thus  P  can  be  expressed  as 

P  =  I2rL.  (8) 

In  addition,  according  to  the  law  of  energy  conservation,  P  can 
also  be  expressed  as 
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P  =  £(Q,h-Q,c).  (9) 

i=  1 

Combining  Eqs.  (5)-(9),  we  get 

(£  rL  =  £  ~  7M  •  a°) 

Combining  Eqs.  (3),  (4),  and  (10),  we  will  get  a  cubic  equation 
group  of  Tii  and  Ti2,  it  is  so  complicated  that  no  analytical  solution 
can  be  found. 

Since  every  TE  module  has  the  same  internal  resistance  r,  in  or¬ 
der  to  get  the  approximation  solution  of  parallel  TEG,  the  following 
approximate  derivation  can  be  done  (see  Appendix  A).  The  parallel 
TEG  which  composed  of  n  TE  modules  can  be  seen  as  n  indepen¬ 
dent  “small  TEGs”,  each  “small  TEG”  works  on  its  respective  heat 
source  and  heat  sink,  the  load  resistance  of  each  “small  TEG”  is 
the  n  times  as  large  as  that  of  the  parallel  TEG,  that  is  nrL,  the 
sum  of  currents  in  “small  TEGs”  is  equal  to  the  current  of  the  par¬ 
allel  TEG.  Therefore,  to  the  ith  TEG,  the  temperature  difference  be¬ 
tween  the  hot-side  and  the  cold-side  of  the  thermocouples  AT,  and 
the  current  I,  are  [24] 

AT,.  =  r,-i-ri2  =  ^E— -  (T,.h  —  t-,0,  (ii) 

r+re+nrL  '  * 


contact  with  the  heat  source  and  heat  sink  by  production  process 
control.  Consequently,  they  have  the  same  thermal  contact  resis¬ 
tance  Rt h  and  Ric ,  according  to  Eqs.  (1)  and  (2),  it  can  be  thought 
that  all  TEG  modules  have  the  same  thermal  conductance  on 
hot-side  and  cold-side.  Therefore,  under  this  ideal  condition,  as 
all  TE  modules  of  parallel  TEG  have  the  same  working  conditions 
and  inherent  parameters,  the  process  and  result  of  performance 
analysis  can  be  simplified. 

Let  Th  be  the  heat  source  temperature,  Tc  be  the  heat  sink  tem¬ 
perature,  I<h  be  the  total  hot-side  thermal  conductance  of  a  TE 
module,  I<c  be  the  total  cold-side  thermal  conductance  of  a  TE  mod¬ 
ule,  and  thus,  Eq.  (16)  can  be  simplified  as 


P  = 


nmoc(Th  -  Tc) 
m2a2B  +  A(r  +  re  +  nrL) 


ru 


(17) 


where 


A  = 


JC  K_ 


(18) 


(19) 


Let  f!^  =  0,  the  maximum  output  power  Pmax  and  match  load  rm 
of  the  parallel  TEG  are  got  as 


nvxATi  _  moc(Tih  -  Tic) 
r  +  re  +  nr  ~  m2oc2Bi  +  Afr  +  re  +  nrL)  ’ 


nm2a2(Th-Tc)2 

(12)  max  4m2a2AB  +  4(r+r^A2’ 


(20) 


where 


At 


K_ 

K^  +  Wc’ 


p  __  T th  Tjc 

'  Kic  Km' 


(13) 

(14) 


r  +  re  m2a2B 

rm=— +t- 

Similarly,  Eq.  (15)  can  be  simplified  as 

m2a2B+A(r+re)  ,  r 
nA  ^  11 


(21) 


(22) 


Therefore,  the  total  current  of  the  parallel  TEG  is 


'  =  £ 


moc{Tih  -  Tjc) 


P  = 


m2u.2Bj  +  A(r  +  re  +  nrL) ' 

The  output  power  of  the  parallel  TEG  is 

m<x(Tih  -  Tic) 


E  m2cc2Bi  +  Ai(r  +  re  +  nrL) 


rL- 


(15) 


(16) 


When  the  model  and  parameters  of  TE  modules  are  selected, 
the  heat  source  temperature  and  the  heat  sink  temperature  of  each 
TE  module  are  known,  and  the  thermal  contact  resistances  are 
measured,  let  =  0,  the  maximum  output  power  and  match  load 
of  the  parallel  TEG  can  be  got.  However,  as  the  Eq.  (16)  is  very  com¬ 
plicated,  it  is  impossible  to  get  the  analytical  expression  of  maxi¬ 
mum  output  and  match  load.  In  practical  application,  when  the 
values  of  all  parameters  are  determined,  the  maximum  output 
power  and  match  load  of  parallel  TEG  can  then  be  calculated  by 
above  method. 


In  Eq.  (22),  rL  is  the  load  resistance  of  the  parallel  TEG, 
^(Th  -  Tc)  can  be  seen  as  the  open-circuit  voltage  of  the  parallel 
TEG,  m2a2g^(r+re)  can  be  seen  as  the  internal  resistance  of  the  paral¬ 
lel  TEG. 

However,  to  the  one  module  TEG,  its  current  is  [24] 


J=  'fiJb-Tc) 
m2oc2B+A(r+re )  ,  ^ 
A  ,L 


(23) 


Thus,  to  the  one  module  TEG,  ^(Th  -  Tc)  is  the  open-circuit 
voltage,  m2a2B*A(r+re)  is  the  internal  resistance,  rL  is  the  load 
resistance. 

Compare  the  Eq.  (22)  with  (23),  it  can  be  seen  that  the  open 
-circuit  voltage  of  the  parallel  TEG  is  as  large  as  that  of  the  one 
module  TEG,  and  the  internal  resistance  of  parallel  TEG  is  n  times 
as  small  as  that  of  the  one  module  TEG.  That  is  to  say,  when  all 
the  modules  have  the  same  parameters,  and  work  in  the  same  con¬ 
dition,  the  parallel  properties  of  the  parallel  TEG  are  the  same  as 
the  parallel  properties  of  common  DC  power,  it  conforms  to  the 
results  of  Zhang  et  al.  [20]. 


3.  Discussion 

3.1.  Output  performances  under  the  ideal  condition 

In  some  cases,  it  is  possible  that  all  TE  modules  of  parallel  TEG 
have  the  same  working  conditions  and  inherent  parameters.  The 
radioisotope  TEG,  which  is  used  in  the  fields  of  military  and  space 
exploration  serves  as  a  good  example,  it  can  be  guaranteed  that 
each  TE  module  has  the  same  heat  source  and  heat  sink  by  design, 
and  it  can  also  be  guaranteed  that  each  TE  module  has  good 


3.2.  Influence  of  contact  effect  on  output  performance 

The  contact  resistance  re  is  the  inherent  parameter  of  TE  mod¬ 
ule.  The  thermal  contact  resistance  has  two  types  of  parameter. 
One  is  the  inherent  parameters  of  TE  module,  I<lt  I< 2,  Ri,  and  R2. 
The  other  is  the  structural  parameters  of  parallel  TEG  which  inflect 
the  contact  conditions  between  the  TE  module  and  the  heat  source 
or  heat  sink,  namely  Rih  and  Ric. 

From  Eqs.  (15)  and  (16),  the  existence  of  contact  resistance  re  is 
just  like  the  increase  of  the  TE  module’s  internal  resistance,  which 
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may  lead  to  the  deceases  of  current  /  and  output  power  P.  There¬ 
fore,  the  effect  of  re  on  output  power  cannot  be  neglected,  though 
its  value  is  generally  very  small  to  commercial  TE  module. 

The  increase  of  thermal  contact  resistance  means  the  decreases 
of  I<i h  and  Kic,  according  to  Eqs.  (13)  and  (14),  and  Bt  will  increase 
with  the  increase  of  thermal  contact  resistance.  According  to  Eq. 
(11),  the  increases  of  A,  and  Bt  will  lead  to  the  decrease  of  AT/,  it 
means  the  decrease  of  temperature  difference  between  the  two 
sides  of  the  thermocouples.  Consequently,  according  to  Eqs.  (15) 
and  (16),  the  current  I  and  output  power  P  will  decrease.  In  a  word, 
the  thermal  contact  resistance  reduces  the  output  power  by 
decreasing  the  temperature  difference  between  the  two  sides  of 
the  thermocouples. 

Overall,  the  existence  of  the  contact  resistance  and  thermal  con¬ 
tact  resistance  will  lead  to  the  decrease  of  parallel  TEG’s  output 
power.  The  result  is  the  same  as  that  of  the  one  TE  module  TEG 
[24]  and  series  TEG  [21]. 

In  extreme  cases,  when  the  contact  resistance  and  thermal 
resistance  of  hot-side  and  cold-side  are  so  small  that  they  can  be 
neglected,  thus,  re  is  zero,  and  I<ih  and  I<ic  tend  to  infinity,  so  Eq. 
(16)  can  be  simplified  as 


P  = 


n 


E 


ma(rih  -  Tic) 

r  +  nrL 


2 


rL. 


(24) 


Let  JE  =  0,  the  maximum  output  power  Pmax  and  the  match  load 
rm  are  got  as 


_m2a2 
max  4  nr 


E(T''h-r>'c 


(25) 


Fig.  4.  Current  I  and  output  power  P  as  a  function  of  load  rL. 


and  heat  sink  structure  is  very  necessary  to  reduce  the  thermal 
resistance  of  heat  sink. 

In  the  practical  application,  improving  the  surface  finish  and 
flatness  of  heat  source  and  heat  sink,  increasing  properly  the  pres¬ 
sure  between  the  TE  module  and  the  heat  source  or  heat  sink  can 
increase  the  real  contact  area.  Besides,  high-conductivity  thermal 
interface  material  (e.g.,  heat  conductive  silicone  grease,  thermally 
conductive  pad,  phase-change  material)  is  often  used  to  ensure 
good  contact  between  surfaces.  These  methods  will  cause  the  de¬ 
crease  of  Ri h,  Rj c  in  Eqs.  (1)  and  (2),  thus,  the  thermal  contact  resis¬ 
tance  will  decrease. 


4.  Experimental  results 


rm=L  (26) 

A  simple  parallel  TEG  composed  of  two  TE  modules  is  used  to 
analyze  the  influence  of  contact  effect  on  output  power.  The  TE 
modules  are  the  TEP1  power  module  of  Xiamen  Taihuaxing  Trad¬ 
ing  Co.  Ltd.,  its  parameters  are 

m  =  126,  a  =  0.000375  VIC1,  K  =  0.476  WK“\  r  =  3.0fl 

It  is  assumed  that  the  two  TE  modules  have  the  same  contact 
resistance  and  thermal  contact  resistance,  set  the  contact  resis¬ 
tance  re  be  0.1  Q,  the  thermal  conductance  of  hot-side  and  cold- 
side  be  all  10  WK_1.  To  one  of  the  TE  module,  Tlh  is  400  K,  Tlc  is 
300  K,  to  another  module,  T2 h  is  450  K,  T2c  is  320  K. 

In  range  of  0-3  Q,  adjust  the  load  resistance  of  parallel  TEG,  and 
the  current  and  output  power  are  calculated  under  two  conditions 
of  neglecting  the  contact  effect  and  considering  the  contact  effect. 
The  result  is  displayed  in  Fig.  4. 

It  is  easy  to  see  from  Fig.  4  that  the  existence  of  contact  effect 
will  reduce  the  current  and  output  power,  therefore,  in  order  to  im¬ 
prove  the  output  power,  the  influence  of  contact  effect  must  be  re¬ 
duced  as  far  as  possible. 

In  the  design  and  production  of  TE  module,  using  high  conduc¬ 
tivity  materials  (usually  copper)  as  the  material  of  conductive  tabs 
is  a  feasible  method  for  reducing  properly  the  resistance  of  conduc¬ 
tive  tabs,  it  can  reduce  the  contact  resistance  of  TE  module.  In  order 
to  reduce  the  thermal  contact  resistance,  a  common  method  is  to 
use  high-conductivity  ceramic  materials  (usually  alumina  ceramic) 
as  the  substrate  of  TE  module,  in  the  condition  that  the  ceramic 
substrate  is  strong  enough,  reducing  properly  its  thickness  will 
make  the  thermal  contact  resistance  decrease.  Besides,  according 
to  the  different  working  conditions,  using  suitable  cooling  method 


In  order  to  verify  the  parallel  model  established  in  this  research, 
a  parallel  TEG  experimental  system  is  built,  the  schematic  of  the 
parallel  TEG  performance  testing  system  is  shown  in  Fig.  5. 

The  parallel  TEG  consists  of  two  TE  modules  which  are  num¬ 
bered  as  TE-1  and  TE-2,  respectively.  The  TE  module  are  the  TEP1 
power  module  of  Xiamen  Taihuaxing  Trading  Co.  Ltd.,  its  inherent 
parameters  are  given  in  Section  3.2  of  this  paper.  The  heat  source  of 
each  TE  module  is  a  heating  plate  made  of  aluminum,  its  temper¬ 
ature  can  be  adjusted  by  the  voltage  on  the  electric  heating  tube. 
The  heat  sink  of  each  TE  module  is  a  water-cooling  plate  made  of 


®  Insulation  material  (2)  Water  cooling  plate  ©  Heating  plate 

©  Electric  heating  tube  ®  Module  TE-1  ©  Module  TE-2 

©  Cooling  water  ®  Voltmeter  (9)  Ammeter 

©  Variable  resistor  ©  E-type  thermocouple 

Fig.  5.  Schematic  of  the  parallel  TEG  performance  testing  system. 
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Fig.  6.  Output  power  P  as  a  function  of  load  rL  under  different  heat  source 
temperature. 


copper,  its  temperature  can  be  adjusted  by  the  flow  and  tempera¬ 
ture  of  cooling  water. 

Some  U-shaped  slots  are  cut  on  the  heating  plates  surface  and 
water-cooling  plates  surface  which  contact  with  the  TE  modules, 
the  dimension  of  slots  is  1.0  mm  x  1.0  mm.  1.0  mm-diameter  E- 
type  thermocouples  (nickel-chromium/copper-nickel)  are  in¬ 
serted  into  these  slots  to  measure  the  temperatures,  the  accuracy 
of  measurement  is  ±0.1  K.  To  reduce  the  thermal  contact  resis¬ 
tance,  the  interfaces  of  the  heat  source  and  the  TE  module,  the 
TE  module  and  heat  sink  are  all  coated  with  heat  conductive  sili¬ 
cone  grease,  it  can  be  thought  that  all  the  interfaces  have  good  con¬ 
tact  addition.  At  the  same  time,  it  is  assumed  that  thermal 
resistance  of  heat  sink  is  so  small  that  it  can  be  neglected,  so  the 
thermal  conductance  I<h  and  I<c  are  all  10  WK_1  which  mainly  con¬ 
sider  the  thermal  resistance  of  ceramic  substrate.  To  minimize  the 
bad  effects  of  convection  and  radiation,  the  TE  module,  the  heat 
source  and  heat  sink  are  all  wrapped  with  insulation  material 
together. 

During  the  test,  the  heat  source  temperature  and  heat  sink  tem¬ 
perature  of  module  TE-1  are  kept  to  400  K  and  300  K,  respectively. 
The  heat  source  temperature  of  module  TE-2  has  the  different  va¬ 
lue,  350  K,  400  K  and  450  K,  while  its  heat  sink  temperature  is  kept 
unchanged  to  300  K.  In  the  different  heat  source  temperature  of 
module  TE-2,  the  value  of  load  rL  is  adjusted  from  0.1  Q  to  5  Q, 
at  the  same  time,  the  output  power  of  the  TEG  system  is  carefully 
measured.  Fig.  6  displays  the  output  power  of  the  parallel  TEG  as  a 
function  of  load  under  different  heat  source  temperature. 

As  shown  in  Fig.  6,  the  results  derived  from  the  model  are  basi¬ 
cally  consistent  with  the  experimental  results.  However,  the  theo¬ 
retical  value  is  slightly  larger  than  the  experimental  results,  the 
deviation  between  the  theory  and  experimental  results  may  be 
due  to  the  previous  assumption  which  ignores  the  heat  conduction 
between  the  thermocouples.  The  temperature  difference  between 
the  two  sides  of  the  thermocouples  will  decrease  when  the  heat 
conduction  between  the  thermocouples  is  considered,  then  results 
the  decreases  of  output  power. 

5.  Conclusion 

A  novel  analytical  model  of  parallel  TEG  was  developed  based 
on  thermodynamic  theory,  semiconductor  thermoelectric  theory, 
and  law  of  conservation  of  energy,  the  equations  of  output  power 
and  current  of  parallel  TEG  were  derived.  The  model  takes  into  ac¬ 
count  much  more  parameters  than  the  one  TE  module  TEG,  such  as 
temperature  of  heat  source  and  heat  sink,  contact  resistance,  and 


thermal  contact  resistance  of  each  TE  module  in  parallel  TEG.  The 
results  show  that  only  all  TE  modules  in  the  parallel  TEG  have 
the  same  inherent  parameters  and  work  in  the  same  condition, 
the  parallel  properties  of  the  parallel  TEG  are  the  same  as  the  par¬ 
allel  properties  of  common  DC  power.  Besides,  contact  effect  will 
reduce  the  output  power  of  parallel  TEG.  The  existence  of  contact 
resistance  is  just  like  the  increase  of  the  TE  module’s  internal  resis¬ 
tance,  which  leads  to  the  deceases  of  output  power.  The  thermal 
contact  resistance  reduces  the  output  power  by  reducing  the  tem¬ 
perature  difference  between  the  two  sides  of  thermocouples.  In  the 
practical  application,  as  the  parameters  of  TE  module  are  fixed,  the 
main  method  of  reducing  the  influence  of  contact  effect  is  reducing 
the  thermal  contact  resistance  between  the  TE  module  and  the 
heat  source  or  heat  sink. 

The  results  derived  from  the  model  are  basically  consistent  with 
the  experimental  results,  it  demonstrates  that  the  analytical  model 
of  parallel  TEG  works  well  for  predicting  the  output  performance. 
However,  the  theoretical  value  is  slight  larger  than  the  experimen¬ 
tal  results  because  of  the  assumption  which  ignores  the  heat  con¬ 
duction  between  the  thermocouples. 

The  analytical  model  of  parallel  TEG  is  a  useful  tool  for  those 
who  want  to  design  parallel  TEG,  since  it  will  give  them  an  idea 
about  how  to  predict  the  output  performance  of  parallel  TEG. 
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Appendix  A.  Approximate  derivation  of  parallel  TEG  circuit 

The  parallel  TEG  circuit  can  be  seen  approximately  as  a  simple 
parallel  circuit  shown  in  Fig.  A.l. 

In  the  simple  parallel  circuit,  there  are  n  power  supplies  in  par¬ 
allel,  and  every  power  supply  has  the  same  internal  resistance  r.  In 
accordance  with  the  Kirchhoff  s  law  and  Ohm’s  law,  we  can  get  the 
following  equations: 


n 


I  —  h  +  h  -T  •  •  •  li  +  •  •  •  4s  In  —  U 

i=l 

(A.1) 

IrL  =  Ui-Iir 

(A.2) 

m  In  r 

1 - 1 

Un 

1 _ 1 

'  ^  r 

1 - 1  1 
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1  1 _ 1 

U2  ; 

' - 1  1 - 

r  . 

1  1 - 1 

(  1 _ 1 

- 1=1 

Ui 
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Fig.  A.l.  Simple  parallel  circuit  diagram. 
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where  I  is  the  current  flowing  through  the  load  rL,  Jz  is  the  ith  branch 
current,  LI,  is  the  electromotive  force  of  the  ith  power  supply. 
Reforming  Eq.  (A.2)  yields 


h  = 


Ui-IrL 


Substituting  Eq.  (A.3)  into  Eq.  (A.l)  yields 

j  Eli^-nfrL 
r 

Reforming  Eq.  (A.4)  yields 

U< 


F 


-  nrL 


Combining  Eq.  (A.l)  and  Eq.  (A.5)  yields 

Ut 


r  +  nrL 


(A.3) 


(A.4) 


(A.5) 


(A.6) 


According  to  Eq.  (A.l),  the  current  I  is  equal  to  the  sum  of  total 
branch  current  /„  that  is  to  say,  in  order  to  get  J,  the  value  of  every  /, 
does  not  matter  as  long  as  the  sum  of  total  //  is  unchanged.  Accord¬ 
ing  to  Eq.  (A.6),  the  sum  of  total  is  equal  to  the  current  I  too, 
therefore,  a  hypothesis  that  replacing  /,  with  can  be  made, 
though  If  is  not  necessarily  equal  to  - that  is 


Ii  = 


Ui 

r  +  nrL 


(A.7) 


Eq.  (A.7)  is  a  formula  to  solve  the  current  of  a  simply  circuit 
based  on  Ohm’s  law,  r  is  the  internal  resistance  of  the  power  sup¬ 
ply,  and  nrL  is  the  load  impendence  of  the  simply  circuit.  Therefore, 
in  order  to  get  the  load  current  I  of  the  parallel  circuit  shown  in 
Fig.  A.l,  the  parallel  circuit  can  be  divided  into  n  simply  circuits, 
and  the  internal  resistance  of  every  simply  circuit  is  r,  the  load  is 
nrL. 
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